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Abstract 

We consider scattering of elastic waves on parallel wedge dislocations in the 
geometric theory of defects or, equivalently, scattering of point particles and 
light rays on cosmic strings. Dislocations are described as torsion singularities 
located on parallel lines, and trajectories of phonons are assumed to be the 
corresponding extremals. Extremals are found for arbitrary distribution of 
the dislocations in the monopole, dipole, and quadrupole approximation and 
the scattering angle is obtained. Examples of continuous distribution of wedge 
and edge dislocations are considered. We have found that for deficit angles 
close to — 2ir a star located behind a cosmic string may have any even number 
of images, 2, 4, 6, . . . The close relationship between dislocations and conformal 
maps is elucidated in detail. 
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1 Introduction 



Geometric description of defects in solids has attracted continuous interest in the 
last years. The central role here is played by torsion. This geometric notion was 
introduced by E. Cartan Jl] (see also 0) and, quite interesting, by analogy with 
elastic media. Later the relation of torsion to dislocations was noticed || || (see 
reviews [5-7]), and disclinations in spin glasses or liquid crystals were found to 
correspond to nontrivial curvature of a SO (3)-connection ||. Elastic media with a 
spin structure (Cosserat media |§) containing arbitrary distribution of dislocations 
and disclinations can be consistently described in the framework of general Riemann- 
Cartan geometry where torsion and curvature are identified with surface densities 
of Burgers and Frank vectors, respectively [I(J (a similar approach was discussed 
in |pl|). Generalizations of gravity models to Riemann-Cartan and general affine 
geometry are reviewed in |12[. It should be stressed that geometric description treats 
single as well as continuously distributed defects in a unified manner. 

The advantage of geometric description of defects in solids is twofold. First, 
in contrast to the ordinary elasticity theory this approach provides an adequate 
language for continuous distribution of defects. Second, a mighty mathematical 
machinery of differential geometry clarifies and simplifies calculations. 

From a mathematical point of view geometric theory of defects in solids and the 
theory of gravity with torsion in the Euclidean formulation are the same models. 
The relation between N point particle solution of three dimensional gravity [13-15] 
and wedge dislocations was first analyzed in [|H| [T7| PL (See also |TS|, [HJ.) In |TU| we 
proposed a three dimensional action for describing defects in solids and showed that 
description of an arbitrary distribution of static parallel dislocations naturally arises 
from an action principle. This solution defines nontrivial geometry on R 3 manifold of 
zero curvature with N removed parallel lines corresponding to torsion singularities. 
Each of the removed lines defines a conical singularity in the perpendicular plane. 
In spite of the fact that the curvature is zero the whole space is not a manifold with 
absolute parallelism because parallel transport of a vector around conical defect 
produces rotation on a deficit angle. A gauge approach to the theory of defects in 
solids and cosmic strings is adopted in ^(J . Recently a series of articles dealing with 
the physical properties of the N point particle solution appeared [21-26]. 

In the present paper we analyze scattering of elastic waves on parallel wedge dis- 
locations. In the eikonal approximation or geometric optics |^7j elastic deformations 
described by the wave equation move along extremals. By analogy with motion of 
photons in electrodynamics it is natural to assume that extremals are trajectories 
of phonons in elastic media with defects Q. In this way analysis of extremals yields 
a complete picture of the scattering of phonons on dislocations. In the case under 

1 There wedge dislocations are called disclinations. The last name is also used in describing 
defects in a spin structure of a media. The geometric theory describes both types of defects, and 
therefore we reserve the name disclination for a defect in a spin structure. 

2 An alternative point of view is adopted in |28| where a point particle is assumed to move along 
geodesic (there extremals and geodesies are called geodesies and autoparallels, respectively). 
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consideration the scattering problem is not the standard one because in general the 
space is not asymptotically flat. For example, the notion of a falling beam of par- 
allel moving particles depends on the distance to the dislocation, and the standard 
definition of a cross section does not work. At the same time one can easily define 
the scattering angle considering asymptotics of a particle trajectory. For one wedge 
dislocation we derive the following scattering angle 

x = — d , (i) 

where 2n8 is the deficit angle that is the angle of the wedge which is removed 
(-1 < 9 < 0) or added (0 < 6) to the media. 

Extremals for one edge dislocation were analyzed perturbatively in |24| . Recently 



extremals were analysed for one dispiration defect containing a wedge dislocation as a 
special case |[25|| . In the present paper we extend these results and analyze extremals 



nonperturbatively for arbitrary distribution of wedge parallel dislocations. In fact, 
a general explicit expression for an extremal is obtained. Already for one wedge 
dislocation we obtain a very rich and interesting scattering picture of phonons. In 
particular, if the deficit angle is close to — 2-k or 9 « —1 then a phonon makes 
several turns around the dislocation before going to infinity. If dislocations are more 
concentrated in a compact domain of the perpendicular plain than at large distances 
one may use perturbation theory. In this way we analyzed monopole, dipole, and 
quadrupole contributions. The dipole approximation coincides with the scattering 
on an edge dislocation. 



The spatial part of the metric for a static cosmic string |29|, |3(J coincides with 
the metric for a wedge dislocation although it was obtained in a quite different 
manner by solving four dimensional Einstein equations for a cosmic string in the 
linear approximation. Therefore scattering of elastic waves on dislocations is the 
same as scattering of point particles or photons on cosmic strings. We derive a 
general formula for the angular separation between images of a star located behind 
a cosmic string. For small angles the expression coincides with the result obtained 
earlier, and there are only two images. If the deficit angle approaches — 2ir then the 
number of images increases because an observer sees the light rays which make no 
rotation around the string, one rotation, two rotations, and so on. 

The results are applicable to a radial disgyration defect in 3 He-A, described by 
a conical singularity, too |JT] . 



In Section 0, a general framework for the analysis of extremals is briefly de- 
scribed. It is applied to extremals in the case of one wedge dislocation in Section 
|3|. In Section f| we calculate the angular separation between the images of a star 
located behind a cosmic string. The analysis of extremals is extended to arbitrary 
distribution of wedge dislocations in Section ||. In Section ^| we show that parallel 
linear dislocations are described by the conformal map of the Christoffel-Schwarz 
type. The dipole approximation corresponding to pure edge dislocation is considered 
in detail in Section |7|. In Section || we study the quadrupole approximation. Two 
simple examples of continuous distribution of dislocations are considered in Section 
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2 A general framework 



We suppose that elastic media with dislocations is an R 3 manifold with some punc- 
tures and/or removed lines corresponding to the cores of dislocations and a given 
Riemann-Cartan geometry [ IT3| . An arbitrary curvilinear coordinate system is de- 



noted by x M , n = 1,2,3. A Riemann-Cartan geometry is defined by metric and 
torsion or, equivalently, by triad and 5'0(3)-connection. Here we use the more con- 
venient second set of variables called Cartan variables. Pure dislocations correspond 
to zero curvature 

R^ ij = dpujl - - (ji <-> v) = 0, (2) 

but nontrivial torsion. Here cJ M y = —ui^ 1 , . . = 1,2,3, is an 5'0(3)-connection 
considered as an independent variable. For zero curvature it is a pure gauge and can 
be always set to zero at least locally by means of local SO (3) rotation. This can be 
done in domains with the trivial fundamental group which are spaces with absolute 
parallelism. If the space contains singularities and one has to remove infinite or 
closed lines from R 3 then this property is violated but may be restored by making 
appropriate cuts. A triad field defines a metric in the usual way 

SV = e^eJSij, Sij = diag(+ + +). (3) 

In our case torsion tensor, 



- ■ 



Tpv = d^ej - u^ 3 e V j -(fjL<-+v), 

is entirely defined by the triad because the SO (3)-connection is set to zero. 

Of course, for a given metric @ one can construct Christoffel's symbols and the 
corresponding curvature as if torsion equals zero. It will be not the same curvature 
as given by @ and is nontrivial in the presence of dislocations. 

We assume that phonons in the media with dislocations move along extremals 
defined by the metric (|3]). Extremals satisfy the well known set of equations 

x» = -T v ^x u x p , (4) 

where T vp p denotes Christoffel's symbols (the Levi-Civita connection) constructed 
from the metric (0). Note that extremals coincide with geodesies in Riemann- 
Cartan geometry if and only if torsion is zero or antisymmetric in all three indices. 
In three dimensions the canonical parameter t is interpreted as the time. Indeed, 
extending three-dimensional space to four- dimensional space-time, x^ 1 — > 
and assuming the following form of the four-dimensional metric 




one finds the fourth equation for an extremal 

x° = x° = at + b, a,b — const. 
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Thus up to a linear transformation the canonical parameter coincides with the real 
time x°. 

Analysis of extremals becomes simpler if one notices that there is always one 
integral to equations (f|) 

g lLV x^x v = C = const > 0, (5) 

which means that the length of the tangent vector or the square of the velocity of a 
phonon (the kinetic energy) remains constant. Another simplification occurs if the 
metric admits a Killing vector field k 11 , 

V^K + V„fc„ = 0, 

where V M denotes the covariant derivative with respect to Christoffel's symbols. In 
this case there exists another integral 

g^k^x" — C\ — const, (6) 

which is related to the symmetry of the space. In two dimensions the existence of 
one Killing vector field is enough for the complete analysis of extremals because one 
has two integrals (f|), (||) for two unknown functions. In higher dimensions complete 
analysis requires the existence of more Killing directions or some other techniques. 
The latter is also required in two dimensions in the absence of a Killing vector field. 

Another approach to the analysis of extremals is provided by the action principle. 
It is well known that equations (f|) follow from the Lagrangian 

L = \g^x\ (7) 

which equals the kinetic energy of a point particle of unit mass moving in a space 
with nontrivial geometry. Of course, here the canonical parameter and real time are 
naturally identified. We see that both phonons and point particles move along the 
same trajectories. The difference is only in the velocity of movement. If a particle 
carries charge or spin which directly interacts with the media then one has to add 
a potential term to the Lagrangian (|7p. In that case a trajectory will differ from an 
extremal. Although the potential term in (0) is absent the motion of a particle is 
nontrivial because the metric explicitly depends on a position. Then equations (f|) 
are nothing more than the Newton law with the force quadratic in the velocity. Due 
to this circumstance the motion of phonons in the presence of a dislocation differs 
drastically from the Newtonian motion of a point particle in a potential field. For 
example, we shall find later a new type of closed trajectories. 

If there is not enough symmetry for the complete analysis of extremals then 
the Hamilton-Jacobi equation may be useful. The Lagrangian (^) gives rise to the 
Hamiltonian 
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where g^ u is the inverse metric and p^= dL/dx fl is the canonical momenta. Then 
the Hamilton-Jacobi equation for the characteristic function Six 11 , t) takes the form 



dt^2 9 dx»dx» ' { } 

or 

Six 11 , t) = W(x^) -Et, E = const, 

where 

1 dW dW 

V" o r = e - (9) 

2 y dx» dx v y ' 

Here E denotes the total (in our case kinetic) energy of a phonon. This equation is 
solved in Section |^ for arbitrary distribution of iV wedge dislocations characterized 
by one Killing direction. 



3 Extremals for a wedge dislocation 

Let us analyze extremals in the simplest case of the media with one wedge disloca- 
tion, characterized by the deficit angle 9 normalized on 27r, so that the real angle is 
equal to 27r# (see Fig.|I]). In this section we derive the formula (0) for the scattering 
angle. The x 3 axis is chosen along the core of dislocation. The corresponding metric 



2ti|I 



Figure 1: A wedge dislocation with the deficit angle 2tt6. For negative 6, — 1 < 6 < 0, 
the wedge is removed from the media. 

has the form |I(J 

ds 2 = dl 2 + (dx 3 ) 2 , (10) 
where the nontrivial two-dimensional part of the metric in polar coordinates reads 

dl 2 = r 2 \dr 2 + r 2 dip 2 ). (11) 

This is the well known form of the metric for a conical singularity with the deficit 
angle 2ir8. It is not asymptotically flat, and particles at infinity cannot be considered 
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as free. They feel a wedge dislocation even at infinite distance. Christoffel's symbols 
corresponding to the metric (|llD are easily calculated, 



-p f — p <p 

rip 1 <pr 





~ ■ 

r 



-(l + 0)r, 
1 + 9 



the others being identically zero. Then equations for an extremal (§) reduce to 

Q 

f = — r 2 + (1 + 0)r(p 2 , 
r 

(p = —2 np, 



x 



0. 



(12) 

(13) 
(14) 



The last equation is the consequence of the translational symmetry along the x 3 
axis, and shows that phonons move along it with a constant velocity. Therefore the 
scattering reduces to the two-dimensional problem in the r, (p plane. The metric 
for a wedge dislocation possesses a rotational symmetry around the x 3 axis. The 
corresponding Killing vector field in polar coordinates has a simple form 




Thus integrals (l|) and (W) for a wedge dislocation reduce to 



r 29^2 + r 2(l+0)^2 



Cn>0, 



r 2(l+^ 



(15) 
(16) 



They have direct physical interpretation to be discussed later. For nonradial motion 
one easily finds the general form of an extremal 



2(1+0) 



where 



sin 2 [(l + 0)(^ + ^ o )] = l, 



(po = const. 




(17) 



18) 



We see that the form of an extremal is parametrized by two arbitrary constants. 
The constant r m is positive and defines the length scale, while the constant <p cor- 
responds to rotational symmetry around the core of dislocation and may be arbitrary. 
Formula (0) describes an infinite number of disconnected branches depending on 
the range of the angle. One branch corresponds to the variation of the sin argument 
from nix to (n + l)ir, |n| = 0, 1, . . . The radius is bounded from below r m < r < oo, 
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and every branch starts and ends at infinity. For n = the angle <p + (fo varies 
from to 7r/(l + 9). For positive 9 the variation of the angle is less than tt, and a 
phonon is repelled from the dislocation. For negative 9 phonons are attracted, and 
can make several turns around the dislocation. Suppose the extremal makes exactly 
m rotations. Then the angle ip + ipo varies from to 27rm, and the corresponding 
deficit angle 9 m is defined by the equation 

(9 m + l)27rm = tt. 

Extremals for a wedge dislocation are drawn in Figs.|2|-|5] for different values of 
9 and <po — 7r 0. In the case 9 > the extra wedge of media of angle 2n9 is added 
to the media. For positive 9 every extremal has two asymptotes going through the 




Figure 2: Extremals for a wedge dislocation with 9 > (i) and — 1/2 < 9 < (ii). 
In (i) we draw two extremals and their reflections for the same 9 but different r m . 
In (ii) we draw two extremals for the same r m but different 9. For 9 = —1/2 an 
extremal makes one turn around the dislocation. 

center and intersecting at the scattering angle x ( see Fig-@(%0- The scattering angle 
is equal to tt — A(p where Atp is the angle between two asymptotes for an extremal 
defined by equation (|17|), (1 + 9)Aip = tt. In this way one immediately gets the 
expression 

7T Tl9 

x 1+9 1+9 

which coincides with (]I|). The scattering angle has clear geometric interpretation. 
It equals one half of the deficit angle times the compression coefficient 

2n9 1 



3 All extremals in this and the following sections are drawn numerically, and the scale is shown 
in the drawings. 
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It means that the wedge dislocation is made by addition of uncompressed media and 
afterwards the media compresses. If one added the compressed wedge of the media 
then the compression coefficient would be absent. If 9 goes to infinity then any 
extremal returns back from the dislocation along the same line. If the dislocation 
is absent, 9 = 0, then there is no scattering and x — 0- It is quite interesting that 
for positive 9 there are points which cannot be joined by an extremal. For example, 
a point on the negative half of the x-axis cannot be joined by an extremal with 
arbitrary point lying in the wedge of inserted media (the shaded region in Fig.^f i J) . 
The reason for this is the conical singularity at the origin of coordinate system. 

For negative 9 the wedge is removed from the media. In this case extremals have 
no asymptotes, but the scattering angle is still defined by the equation ([!]). In Fig. || 
(ii) we show two extremals for different 9 but the same r m = 1 which have up to 
one turn around the dislocation, one turn corresponding to 9 = —1/2 when half 
of the media is removed. If— 1 < 9 < —1/2 then the phonon makes several turns 
around the dislocation before going to infinity as shown in Figs. |3|-[5[ This kind 




6=-2/3 
X=-2n 



Figure 3: An extremal for 9 = —2/3 makes one turn around the dislocation and 
then goes forward in the same direction. 

of trajectories is not specific only to the dislocations. Similar extremals are met 
already in the Schwarzschild space-time |32| . These trajectories produce multiple 
images of a star located behind the singularity. 
At infinity r — > oo equations ([15|), fllCf ) yield 

-^—r 1+e = jQjt + const. (19) 
1 + 9 v v ' 

It means that for 9 > — 1, a particle moves to space infinity in infinitely long time, 
and the manifold is complete at r — > oo (any extremal can be continued to infinite 
value of the canonical parameter in both directions). Later we show that the radial 
extremals can be continued through the singularity and become complete too. In 
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Figure 5: An extremal for 9 = —4/5 makes two turns and a half around the wedge 
dislocation before going to infinity. 



this sense the whole manifold is complete. Unphysical region 6 < — 1 corresponding 
to the case when one removes more than the whole media leads to an incomplete 
manifold which should be extended. 

Physical interpretation of the first integrals (ji5|), ( |16D is clear from the following 
consideration. In the primed coordinate system, 



1 



.1+0 



1 + 9 

the line element (|Tl|) becomes flat 



dt z = dr' 2 + r ri dip 



J'2 



This means that in these coordinates phonons move along straight lines, and the 
energy, momentum, and angular momentum are conserved. The kinetic energy of a 
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phonon of unit mass is easily calculated 

E'= 1 -(r- + r^)= 1 -C . 

So the first integral fll5|) is proportional to the kinetic energy or the square of the 
velocity. The momentum is defined by Cq and the direction of the trajectory. The 
calculation of the angular momentum, 

yields physical interpretation for the other integration constant. It is defined by 
the initial condition of the system and is arbitrary. The constant of motion Cq 



exists for the wave equation of acoustic waves [27] and is different for longitudinal 
and transverse phonons. This is the constant distinguishing trajectories of different 
waves. One may assume that point particles also move along extremals if they do 
not carry charge or spin which directly interact with the media. Then the constant 
Cq is parametrizing the type of a phonon and a point particle. For a point particle 
it is defined by the initial condition and may be arbitrary. The form of an extremal 
( |TTD is defined by the constant (|18|) depending both on Cq and C\. The constant 



C\ is arbitrary and therefore phonons and point particles move along the same 
trajectories, the only difference being the velocity of movement. 

Suppose that the external observer examines the scattering using some device 
which does not interact with the media. Then he observes the scattering on dislo- 
cation in polar coordinates r, (p which are flat for his device. For him the existence 
of dislocation exhibits itself as the bending of phonon trajectories. 

Let us compute physical characteristics of the scattering phonon as measured 
by an external observer for whom a dislocation appears through the properties of a 
phonon. The observed kinetic energy, 

E= 1 -(r 2 + r^) = 1 -CQT-^, (20) 

is positive and explicitly depends on the distance to the dislocation. It decreases 
and increases for positive and negative deficit angle 9, respectively, as the phonon 
approaches the dislocation. At the infinity the energy tends to zero for positive 9 
and to infinity for negative 9: 

9>0: E 0, r ^ oo, 
9 < : E — > oo, r ^ oo. 

Of course in a real world the energy cannot increase to infinity and is bounded by 
the finite size of the bodies. The radial and angular components on the momentum 
have the form 



p r = f = Jc~ o r- 2e ~ 1 Vr^+ () )-r 2 Ji 1+0 \ (21) 
Pip = r^ = C 1 r- 2e - 1 . (22) 
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The absolute value of the angular component increases and decreases as the phonon 
goes closer to the dislocation for 9 > —1/2 and — 1 < 9 < —1/2, respectively. 
For 9 = —1/2 it does not depend on the distance. The behavior of the absolute 
value of the radial component is more complicated. For — 1 < 9 < it grows up 
monotonically from zero to infinity as the radius grows from r m to oo. For 9 > it 
is zero both at r m and infinity and has one maximum at 



l + 2#\ 2(1+9) 



29 J 

The angular momentum as measured by an external observer, 

L = r 2 Lp = C ir - 2e , 

also depends on the distance. 

Radial extremals must be treated separately. Equation (|I6|) for <p = const and 
Ci = is satisfied, and equation (jl6|) can be easily integrated 



r l+e = y/Co(l + 9)(t + t ), t = const. 

For any deficit angle the core of dislocation is reached by the radial extremal at 
finite time. It is naturally continued through the conical singularity: we consider two 
radial extremals with angles ip and ip+tr as two halves of one complete extremal. This 
is one possibility when the phonon is considered to go through the conical singularity. 
The other possibility is to assume that the phonon reaches the singularity in finite 
time and disappears. The present model does not allow to choose between the 
alternatives, and one has to invoke other theoretical or experimental arguments. 
This is beyond the scope of the paper. 

Note that circle extremals, r = const, are absent as the consequence of equation 
(|P2"D, though integrals (|15D, (p~6|) admit this solution. This happens because to obtain 
the first integrals equations (|12"D, ( |13D must be multiplied by f. 

The analysis of extremals shows that the motion of phonons is not a potential 
one. That is, there is no effective potential yielding the same trajectories in the flat 
Euclidean space. Indeed, the observed kinetic energy ( PD| ) suggests the potential 

U = --Cor~ 2e + const 
2 

for the total energy to be conserved. Then one gets the Lagrangian 

T * 2 i 1 2-2 i 1 ri -20 

L = -r +-r + -C r 

with Euclidean metric but nontrivial potential. This Lagrangian yields qualitatively 
different trajectories for a particle. 
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4 Cosmic strings and multiple images 



The metric for a wedge dislocation with a negative deficit angle, 9 < 0, coincides with 
the spatial part of the metric for a cosmic string first found in |33|] . It is interesting 
that the metric for a cosmic string was obtained in a quite different way by solving 
four-dimensional Einstein equations in the linear approximation. Nevertheless the 
metric is essentially the same, and therefore extremals found in the previous section 
describe also light rays in the presence of a cosmic string. The well known effect of a 
cosmic string is the appearance of double images of a star located behind the cosmic 
string, see Fig. Wi). In this section we derive a general formula for the angular 




Figure 6: Double images of a star located behind a cosmic string with small deficit 
angle (i). If a deficit angle is close to — 2n then a star may have any even number 
of images. In Fig. (ii) we show a star with four images, two of them being produced 
by light rays making one turn around a cosmic string. 

separation between the images. Originally, the angular separation was obtained 
for small angles by calculating the angle in the coordinate system where extremals 
are straight lines (that is, in the media with removed wedge but not deformed). 
We obtain the formula which is valid for all angles and in the coordinate system 
describing the media after the dislocation is created. 

Suppose that a star and an observer are located at distances I and d, respectively, 
as shown in Fig. ^. The angle ip between the tangent vector to a light ray and the 
radial direction is defined by the expression 



tan?/> = 




Using equation (|17|) one easily finds 

V = -(l + W + 0o), 
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up to an addition of 2nn, \n\ = 1,2,... Location points of a star and observer are 
defined by the equation 



2(1+61) 



— sin ijji 



2(1+9) 



sin if) 2 



(23) 
(24) 



where 



ipi = -(1 +6)(tt + 2m7c + ip ), 
i> 2 = -(! + %,„. 



(25) 
(26) 



Here m is the number of rotations of the light ray around the cosmic string before 
reaching the observer. The whole number of rotations of the corresponding extremal 
may be greater. The difference between ip 2 and ipi must be less then n. Otherwise a 
star and an observer would be crossed by different branches of an extremal. There- 
fore the inequality 

(1 + 0)tt(1 + 2m) < 7T (27) 

restricts the maximal number of rotations for a given deficit angle. As the conse- 
quence of equations (ESI) and (Effl) we have 



V>2 = ^i + 7r(l + 0)(l + 2m). 

Denoting the angular separation between images by A = 2{ir— ip 2 ) one easily deduces 
from ( |23"D , (p3|) a general formula for the angular separation between two symmetric 
images produced by the light rays making m rotations around the string 



A 



sin • 



sm 




7r(l + 0)(l + 2m) 



(28) 



Here m is any natural number satisfying the inequality (|27J). For small deficit angle, 
\0\ <C 1, an extremal cannot make a turn around the dislocation, m — 0. Therefore 
there are only two images separated by a small angle, A 1. In this case equation 
( |28|) reduces to the well known formula |29[ 



A 



d + V 



I27T0I. 



(29) 



The coincidence of the results is not trivial, because the angle is computed in different 
coordinate systems. It means that the coordinate transformation describing a defect 
creation must be described by a conformal transformation which preserves angles. 
This is indeed the case as will be shown in the next sections. In fact, this assumption 
was made but not written in the original derivation of (|29|) . 



14 



Let us consider a general case in more detail. Light rays emitted by a star 
are parametrized by the angle ipi. They cross the x axis at point d defined by the 
equation 

d\ 2{1+e) sin 2 [^ + (l + g)7r(l + 2m)] = x 
I J sin 2 ?/'i 

We see that for a given m the observer at any distance d sees two images: one above 
and one below the cosmic string as shown in Fig. [f| If the deficit angle is sufficiently 
close to —2tt then an extremal can make several turns around the cosmic string. 
If the maximal number of rotations for which the inequality ( p7|) is satisfied equals 
M then the observer sees 2(1 + M) images: two images caused by the light rays 
which do not make a rotation before reaching the observer, two images because of 
the extremals making one rotation, and so on up to M rotations. In Fig. ^(ii) we 
show a star with four images for 9 = —3/4. For this deficit angle any extremal 
makes one full rotation around the string but two of them reach the observer before 
making a turn as shown in the picture. Light making an additional rotation reaches 
an observer later and produces a time delay. 



5 Extremals for N wedge dislocations 

Let us consider an arbitrary number N of parallel wedge dislocations, characterized 
by deficit angles 8 n , n = 1, . . . , N. The corresponding two-dimensional metric dl 2 in 
complex coordinates z = x + iy has the form |T0| 



dl 2 = dzdzl[[(z-z n )(z-z n )} e ". (31) 

n 

Here a bar denotes complex conjugation. Dislocations intersect the x,y plane in 
the points z n . Direct integration of the corresponding equations for extremals is 
sophisticated in this coordinate system because in general there is no Killing vector 
field in the x, y plane. Even in the case of one edge dislocation (the dipole of 
two wedge dislocations) the problem requires numerical analysis if the metric is 
first expanded at large distances and then extremals are found |24|]. However, the 
problem becomes almost trivial if one does not expand the metric and note that there 
is a coordinate system where the metric becomes Euclidean. In complex coordinates 
denoted by w it has the Euclidean form 

dl 2 = dwdw, (32) 

where 



w = u + iv = e ia 



[ Z d(l[((-z n ) 9 "+C. (33) 



Here a and C are arbitrary real and complex numbers. The contour of integration 
should not cross and contain closed loops around points z n . We discuss this integral 
in detail in Section |] in several simple cases. If all dislocations go through the x- 
axis, lmz n = 0, then this function is called the Christoffel-Schwarz integral (see, 
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for example, |3"3fl ) and is well known in complex analysis. It conformally maps the 
upper half plane z > onto the polygon with one vertex lying at infinity and other 
vertices at points w n . Each point z n is mapped into the vertex w n with the inner 
angle 7r(l + # ra ). In a general case when lmz n ^ for some dislocations the conformal 
map becomes more involved. 

The integral (133]) yields the solution to the Hamilton-Jacobi equation (R|) 



S = V2E\J (w — c)(w — c) — Et, c = a + ib = const. 
For the Euclidean metric (^) extremals are straight lines, 

{u — a) sin 7 = {y — b) cos 7, (34) 

or 

(w - c)e" i7 = (w - c)e h (35) 

going through a point (a, b) at an angle 7. Scattering of a phonon in real crystal 
occurs in the original Cartesian coordinates Therefore equation (|34] ) has to be 
written in the original system. This cannot be done for arbitrary distribution of 
defects because the integral ( ]33| ) cannot be taken in a general case. To describe 
scattering of phonons one needs only asymptotics at large distances \z\ 3> \z n \. In 
the following sections we analyze the scattering up to the quadrupole expansion. It 
is characterized by three global parameters of a distribution. 
The total deficit angle (or the charge) normalized on 27r is 

N 

e="£6 n . (36) 

n=l 

We assume that > — 1 for the physical reason because one cannot withdraw more 
than the whole media. The total Burgers vector (the dipole momentum) and the 
quadrupole momentum normalized on 2n are 

N 

B = J2 6 nZn, (37) 

71=1 

1 N 11 

M = - E tfnOk ~ OJnk) Z n Z k = ^B 2 --J2 0n4 (38) 
Z n,k=l Z Z n 

In this definition the Burgers vector and the quadrupole momentum should be con- 
sidered as complex numbers. In the next section we define the Burgers vector as a 
vector. It has the same absolute value as B but differs from the vector joining the 
origin of the coordinate system with the point B in the complex plane. 

The integral (|33J) has different asymptotics depending on the total deficit angle. 
We write down the first three terms for different values of B. 

w « e ia (o^y^ e+1 - ^z e + Q^-y^ 1 ) + C, 0^0,1, (39) 
w w e ia (z-Bhiz- —\ +C, 6 = 0, (40) 

w w e ia Qz 2 -Bz + Mlnz) +C, = 1. (41) 



16 



In the case 6 = no wedge of media is added or removed. The values 9 = 1,2,... 
determine the position of ln-term in the expansion of w(z). The quadrupole term is 
logarithmic for = 1. 

Expansions (|39|)-(f4T|) show that any extremal going to infinity in the z-plane goes 
also to infinity in the w-plane. In the latter case extremals are obviously complete, 
and thus they are complete in the z-plane. If an extremal does not go to infinity 
then it is also complete because it is either closed or it can be continued to infinite 
value of the canonical parameter within a finite domain. Note that we have only 
conical singularities and extremals can be naturally continued through them. 

If the total angle 8^0 then the behavior at large distances is characterized by 
the first term in equation (|39|) describing single wedge dislocation discussed in the 
previous section. Therefore we consider the dipole and quadrupole approximation 
in the following sections. 



6 Conformal maps and dislocations 

The integral ( |33"D provides a conformal map between a plane with arbitrary number 
of conical singularities corresponding to wedge dislocations (z-plane) and a Rieman- 
nian surface (covered by u>-coordinates) with the Euclidean metric. This map has a 
clear physical interpretation corresponding to a defect creation. One starts with the 
Euclidean plane, then cuts or adds parts of the plane making a Riemannian surface 
with boundaries. Then the boundaries are glued together. The gluing process is 
uniquely defined by the conformal map. The points on a boundary are considered 
as the same if they correspond to the same point on the cut in the z-plane. Below 
we consider several examples of conformal maps defining the defects considered in 
the present paper. 

The Christoffel-Schwarz integral depends on two constants. We suppose that 
the point z\ always coincides with the origin of the coordinate system in the z-plain, 
and no point z n , n = 2, 3, . . . , lies on the negative half of the real axis. Then the 
constants may be chosen in such a way that the point w± coincides with the origin in 
the u>-plain and negative part of x-axis is mapped onto negative part of the w-axis. 

The simplest example of the Christoffel-Schwarz map (|33p is provided by the 
wedge dislocation discussed in Section |[ In this case we have only one dislocation 
core located in the origin of the coordinate system 

Zl = o, B x = 9. 

The corresponding conformal map is given by the function 

w = —/ + \ (42) 

and is shown in Fig. |7|. The Christoffel-Schwarz integral maps the upper half plane 
into the polygon with two vertices. One vertex is located at infinity while the other 
coincides with the origin w\ = and has inner angle 7r(l + 9). The whole z-plane 
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Figure 7: The Christoffel-Schwarz conformal map for a wedge dislocation. 

with a cut along the positive half of the real axis is mapped onto the w-plane with 
the wedge of deficit angle 2n9 removed or added to the plane for — 1 < 9 < or 
< 9, respectively. This conformal map reflects the process of the creation of a 
wedge dislocation. One takes an infinite elastic media without a defect, cuts out the 
wedge 2ti\9\ for — 1 < 9 < 0, and identifies the points corresponding to the same 
point on the cut in the 2-plane. For < 9 one has to cut the w-plane, move apart 
the sides, and insert the wedge of a media without stresses. 

The dipole of two wedge dislocations is called an edge dislocation and is char- 
acterized by a constant Burgers vector B to be defined later. To construct the 
conformal map we choose the following locations and angles of the wedge disloca- 
tions 

Zl = h a 1 = 7°' ( 43 ) 

z 2 = h, 9 2 = 9. 

We assume that < 9 < 1. Then the Christoffel-Schwarz integral (^) takes the 
form 

w=[ Z d((^). (44) 



c 

Here the constants are chosen in such a way that Wi(zi) = 0. Location of the second 
vertex in the upper half plain (see Fig. ||) is expressed in terms of the gamma- 
functions 

W2 = £d(^^\ = he™ 6 dx (^—^j =he ine T(l-9)T(l + 9). (45) 

In Fig. [5] we show the conformal map for 1/2 < 9 < 1. The z-plane has a cut along 
the positive part of the real axis. The upper half plane is mapped onto the triangle 
in the upper half of the u>-plain with one vertex lying at infinity. The lower half 
of the z-plain is mapped symmetrically. The whole z-plain with the cut is mapped 
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Figure 8: The Christoffel-Schwarz conformal map for a dipole of two wedge dislo- 
cations representing one edge dislocation. 

onto the u>-plain with the removed strip as shown in the picture. We define the 
Burgers vector as the vector joining two points u>_ and w + corresponding to lower 
and upper shores of the cut in the z-plane. For an edge dislocation for x > z 2 it is 
a constant vector. The absolute value of the normalized vector is 

27r|B| = w + — u>_ = 2lmw2- (46) 

This definition is applicable for arbitrary distribution of wedge dislocations. In 
a general case the cut may be defined as the line joining the sequence of points 
Zx, . . . , zn, oo including the infinite point. Of course, the corresponding Burgers 
vector may not be constant. Using the expression ( fl5|) and the property of gamma- 
functions we get for an edge dislocation 

2it\B\ = 2hT(l - 0)r(l + 9) sm(n6) = 2hn6. 

This remarkable result shows that the normalized Burgers vector equals the product 
h9. Note that this is an exact result for a dipole of two wedge dislocations, and it 
coincides with the definition ( j3?|) for the distribution (|E3|) used in the large distance 
expansion of the Christoffel-Schwarz integral. 

Let us compare the conformal map for a dipole of wedge dislocations with the 
conformal map (f40D of the pure edge dislocation 9 = 0, and 5^0 

w = z - Blnz + iirB, (47) 

where we have specified the values of the constants. Elementary analysis yields 
the conformal map shown in Fig. ^|. It maps the upper half plane of the z-plain 
onto the upper half plain of the u>-plain with the cut as shown in the picture. The 
origin of the coordinate system is mapped into the infinite point, wi(z\) = oo. This 
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Figure 9: The conformal map in the dipole approximation (an edge dislocation). 



conformal map has less clear physical interpretation but is more universal because 
it does not depend on the details of the distribution at small distances. 

The pure quadrupole distribution of wedge dislocations may be organized in the 
following way 



zi = 0, 


0i 


= -0, 


z 2 = h, 


02 


= 0, 


Z3 = I, 


3 


= 0, 


= l + h, 


04 


= -0, 



z 4 

where h, I, and 9 are constants. For defmiteness we assume that I > h and < < 1. 
Then only the quadrupole momentum differs from zero 

= 0, B = 0, M = 9lh. 

The corresponding conformal map is provided by the Christoffel-Schwarz integral 

'(C-MC-0 " 



w 



(48) 



It cannot be taken in elementary functions but may be analyzed qualitatively. The 
location of the vertices in the upper half of the w-plain is given by the convergent 
integrals 



w 2 



w 3 



U>4 



he 



in8 



dx 



w 2 + (l- h) 
w 3 + he~ md 



[1 — x){p — x) 
x(p + 1 
1 

dx 



X 



X 



dx 



(q + x)(q + l—x) 
(1 — x)(p — x) 



x(p + 1 



X) 
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where 



We see that I111W4 = 0. The lower half plain is mapped symmetrically. The cor- 
responding conformal map is shown in Fig. [H]. This conformal map has evident 




Figure 10: The conformal map for the four wedge dislocations representing the 
quadrupole source. 

physical interpretation. The dislocation is made by withdrawing a hexagonal tube 
of media. 

In the pure quadrupole approximation, © = 0, B — 0, the conformal map for 
arbitrary distribution of defects ( flOf) takes the form 

M , , 

w = z-—, (49) 

where the constants have been specified. It maps circles into ellipses as shown in 



Fig. [11]. The circles of radius greater than v M are mapped preserving orientation, 
and the outer part of the z-plain, \z\ > \f~M is mapped onto the whole w-plain 
with the cut along the imaginary axis. The inner circles are mapped with opposite 
orientation, and the interior \z\ < also covers the whole w-plain. We see that 
the function (f49|) maps the whole z-plain onto a Riemann surface consisting of two 
copies of w-plain glued together along the cut. We leave it to the reader to imagine 
how this conformal map reflects the defect creation. 



7 The dipole approximation 

If the total deficit angle equals zero, O = 0, and 5^0 then the scattering at large 
distances is described by the dipole term in (|40|). The line element corresponding 
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Figure 11: The conformal map in the quadrupole approximation. 



to the conformal transformation (|47| ) in polar coordinates is asymptotically flat 

2 cos ip N 



dl 2 ={l-B 
The extremals (153) have the form 



(dr 2 + r 2 dtp 2 ). (50) 



r sm(ip — 7) + B(\n r sin 7 + cos 7) — a sin 7 + 6 cos 7 = 0. (51) 

In the dipole approximation the scattering is characterized by the vector, and trajec- 
tories depend on the angle at which extremals go to infinity. We consider two cases. 
For 7 = the trajectories at infinities are perpendicular to the Burgers vector. In 
this case equation (|51|) reduces to 

y . y + b 

— = — tan , 

x B ' 



and extremals are parametrized by one arbitrary constant b. In figure [12] extremals 
are shown for 6 = and b > 0, respectively. Extremals for b < are obtained from 
those depicted in Fig. ^ (ii) by the reflection y — > —y. Extremals have asymptotes 
at infinity shifted by 71 B. So the net result of the scattering is the shifting of a phonon 
trajectory defined by the Burgers vector and a time delay. Note the existence of 
returning trajectories near the dislocation when phonons move from the right. If 
phonons move from the left there is no backward scattering. It means that an edge 
dislocation is "invisible" when it is being seen perpendicular to the Burgers vector 
from the left. 

The extremals which are parallel to the Burgers vector at infinity have different 
behavior and are shown in Fig. [H| They correspond to 7 = 7r/2 and are defined by 
the equation 

x 2 + y 2 = exp (2^) . (52) 
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Figure 12: Extremals for the edge dislocation for b = (i) and 6=1/2 (ii) which 
are perpendicular to the Burgers vector at infinity. For any b ^ the returning 
extremal goes through the core of dislocation. 

In this case they are parametrized by one arbitrary constant a. At infinity the 
extremals have no asymptotes. This means that even at infinity the observer can 
feel the dislocation because trajectories are not straight lines. There is also a time 
delay. For a = a where 

do = B In B 

the extremal has the intersection point at (-5,0). For a > an extremal has one 
branch. If a < clq then solution of equation ( |52| ) has two branches. One branch starts 
and ends at infinity while the other is closed and surrounds the dislocation core. In 
fact this branch may not be closed if the velocity along the x 3 axis is nonzero. In 
that case it looks like a spiral. 

Let us remind the reader that extremals are drawn for the dipole approximation, 
and therefore one may expect qualitative agreement at large distances from the core 
of dislocation. 

Extremals for an edge dislocation were analyzed numerically in ||24|| . Qualitative 
behavior is the same except for the loss of returning trajectories in Fig. [12| and closed 
trajectories in Fig. [13|. 

8 The quadrupole approximation 

In the quadrupole approximation, = 0, B = 0, M ^ 0, the conformal map (|40|) 
at large distances reduces to 

M w 
w = z . (53) 

z 

Without loss of generality we set ImM = and suppose that M > 0. To get the 
scattering for negative M one has to turn the whole picture on the angle ir. The 
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a 2 <a 1 <a 



Figure 13: Extremals for an edge dislocation which are parallel to the Burgers vector 
at infinity. They have no asymptotes. For a < a an extremal has two branches. 

line element corresponding to the conformal map (|53|) has the form 

dl 2 = ( 1 + 2M^^j (dr 2 + r 2 d^ 2 ) (54) 

and goes to the Euclidean line element faster than in the dipole approximation as 
r — > oo. In the quadrupole approximation equation fl3~4"|) for an extremal reduces to 



M 

r sin(<^ — 7) H sin(<£> + 7) — a sin 7 + b cos 7 = 0. (55) 

r 

At infinity r — > 00 this equation coincides with that for a straight line, and therefore 
any extremal has asymptote. Let us consider two cases. Extremals parallel to the 
x axis at infinity correspond to 7 = and are defined by the equation 

x +y = — . (56) 

y + b 

There is no extremal for 6 = 0. For negative b < all extremals lie in the upper 



half plane and are shown in Fig. |14|(%). Extremals for b > are obtained from those 
for b < by reflection y — > —y. For — 1\[M < b < an extremal has only one 
branch going from x = — 00 to x = +00 and touching the x-axis at the origin of 
the coordinate system. If b = —2\fM 1 then the extremal has the intersection point 
at x = 0, y = vM- For b < —2y/M an extremal consists of two branches. One 
infinite branch goes from x = — 00 to x = +00 and does not touch the x-axis. The 
other branch is closed and goes through the origin. Any infinite branch has the same 
asymptote y = — b at both sides as x — > ±00, and the scattering reduces to a time 
delay. 
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Figure 14: Extremals in the quadrupole approximation which are parallel to x— 
(i) and y-axis (%%) at infinity. An extremal has two branches for b < —2\[M and 
arbitrary a. 



The extremals parallel to the y axis at infinity correspond to 7 = tt/2 and are 
defined by the equation 

, n MX 

x 1 +y = . 

x — a 

It differs from (|56]) by the sign in the right hand side and the exchange x «-> y. Here 
M is also supposed to be positive. For a = the extremal is the circle of radius \[M. 
Extremals with positive a > are shown in Fig. Ti(ii). Extremals corresponding to 
negative a < are obtained by the reflection x — > —x. Any extremal consists of two 
branches. One branch is infinite and goes from y = —00 to y = +00 with the same 
asymptote x = a. The other branch is closed and located within the circle on the 
other half of the plane. It goes through the origin. So the scattering reduces only 
to a time delay. 



9 Continuous distribution of defects 

One of the main advantages of geometrical description of defects is the possibility to 
describe continuous distributions of dislocations and disclinations. First, we consider 
circularly symmetric distribution of wedge dislocations. They are assumed to be 
uniformly distributed on the disk of radius R with the following density of deficit 
angles 

f q, r<R, 
10, r> R. 



p(r) = 

The normalized total deficit angle for this distribution is 
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The conformal factor of the corresponding metric 



(57) 



satisfies the three-dimensional Einstein equations which in the considered case re- 
duce to the Poisson equation [TO] 



and is defined by the integral 



The result is 



A0 = —p(r) 



— [ dsp{s) In \s — f\. 
2n J 



qr 



2 4 4 

^-R 2 lnr, 
2 

Outside the dislocations the line element has the form 



r <R, 
r > R. 



(58) 



df 



r« R '(dr 2 + r 2 d(p 2 ). 



This means that the distribution of wedge dislocations as seen from the outside is 
the same as for one wedge dislocation ([□]) with the deficit angle 0. Therefore tra- 
jectories of phonons at large distances are the same as for a single wedge dislocation. 
This is the result one would not expect because Einstein equations are nonlinear, 
and, in general, there is no superposition principle. Note that metric (57) with the 
conformal factor ( |58"D is the exact solution to the Einstein equations. 

Now we consider the continuous distribution of edge dislocations characterized 
by the density of the normalized Burgers vector f3(r)= (/3 x ,/3 y ). It is assumed to be 
directed along the x-axis, (3 y = 0, and uniformly distributed on the disk of radius R 



(3, r<R, 
0, r> R. 



In this case the Einstein equations also reduce to the Poisson equation for the 
conformal factor 

A0 = 2tt(V/3). 



Its solution has the form 



(3r cos <f 



PR 2 cos f 



r < R, 



r > R. 



(59) 
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At large distances, r 3> R, the space is described by the line element ( |5UD with 



B = (3nR 2 , (60) 

and is asymptotically flat. We see that at large distances the continuous distribu- 
tion of edge dislocations behaves like one edge dislocation with the total Burgers 
vector (|60|) . The superposition principle is valid only for the main term in the line 
element. There are further corrections which distinguish one edge dislocation from 
the continuously distributed ones. The first correction may be easily obtained in 
complex coordinates. The conformal factor ( |5UD outside the distribution yields the 
line element 

dl 2 = exp ( —B ( — hi dzdz. 



It may be written in flat form dl 2 = dwdw, where 
w = J d(exp 

This expansion yields the quadrupole correction, M = B 2 /2, to the continuous 
distribution of edge dislocations. 




10 Conclusion 

In the geometrical theory of defects an elastic media is assumed to be a manifold 
equipped with a Riemann-Cartan geometry, curvature and torsion tensors describ- 
ing the surface densities of disclinations and dislocations, respectively. When discli- 
nations are absent then curvature is equal to zero and the corresponding SO (3)- 
connection can be locally gauged away. Then geometry is entirely defined by the 
triad satisfying the three-dimensional Einstein equations. For a given triad (or a 
metric) one has two choices: either to consider a manifold as a Riemannian one 
with Christoffel's symbols as a connection and nontrivial curvature or to introduce 
a torsion in such a way that curvature corresponding to the metrical connection is 
identically zero (a manifold with absolute parallelism if one excludes paths around 
the singularities) . In the last case the manifold has nontrivial torsion defined by the 
triad only. In this sense three-dimensional gravity and geometric theory of disloca- 
tions differ but the solution to one problem yields immediately the solution to the 
other and vice versa. That is a solution to the three-dimensional Einstein equations 
may be always interpreted as some distribution of defects in elastic media. The 
correspondence between the notions is as follows: 

spinless point particle = wedge dislocation 

point particle with spin = screw dislocation 

static point particles = parallel dislocations 

moving point particles = non parallel dislocations. 
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Recently solutions for moving particles with and without spins were found in three- 
dimensional gravity [35-38]. It means that the Euclidean version of this solution 
describes an arbitrary number of arbitrary oriented wedge and screw dislocations. 
This problem is so complicated in the framework of ordinary elasticity theory that 
it can hardly be solved. At the same time it can be solved in the framework of the 
geometric theory of defects. 

Another advantage of the theory of defects concerns the propagation of elastic 
waves in the presence of defects considered in the paper. We are not aware how to 
formulate this problem unambiguously even for a single dislocation in the frame- 
work of elasticity theory. At the same time the geometrical approach adopted in the 
present paper provides an almost standard and simple answer: the perturbations 
propagate along extremals for a metric describing given distribution of dislocations. 
In this way we have found trajectories of phonons. Their behavior is quite different 
from that in the potential motion of point particles. We hope that this non-potential 
motion of phonons in the presence of defects may be confirmed or refuted experi- 
mentally in the future. Maybe elastic media provides a much better experimental 
field for testing gravity models. 

In the present paper we have shown that a star located behind a cosmic string 
may produce any even number of images. The number of images is larger when the 
deficit angle is closer to —2tt. Analogous behavior of extremals is encountered also 
in general relativity already for the Schwarzschild solution [ 32| . This effect leads 
to the appearance of multiple images of a star located behind the singularity and 
changes the usual description of gravitational lensing. If the singularity is strong 
enough then an observer sees not two but a large even number of images of a star. 
There is a wide field of speculations about what we really see in the sky. 

This work is supported by the Russian Foundation for Basic Research, Grants 
RFBR-96-010-0312 and RFBR-96-15-96131. 
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